Phase diagrams of correlated electrons: systematic 
corrections to the mean field theory 

R. Hlubinai'2 

^Division of Solid State Physics, Comenius University, Mlynska Dolina F2, 842 
48 Bratislava, Slovakia 

2 International School for Advanced Studies (SISSA), Via Beirut 2-4, 1-34014 
Trieste, Italy 

Abstract. Perturbative corrections to the mean field theory for particle-hole 
instabilities of interacting electron systems are computed within a scheme which 
is equivalent to the recently developed variational approach to the Kohn-Luttinger 
superconductivity. This enables an unbiased comparison of particle-particle and 
particle-hole instabilities within the same approximation scheme. A spin-rotation 
invariant formulation for the particle-hole instabilities in the triplet channel is 
developed. The method is applied to the phase diagram of the t-t' Hubbard 
model on the square lattice. At the Van Hove density, antiferromagnetic and 
d-wave Pomeranchuk phases are found to be stable close to half filling. However, 
the latter phase is confined to an extremely narrow interval of densities and away 
from the singular filling, d-wave superconducting instability dominates. 



1. Introduction 

In strongly correlated electron systems, usually several symmetry breaking states 
appear as natural candidates for the ground state. This is so because, due to the 
large interaction strength, typically various mean field criteria for an instability are 
simultaneously satisfied. Unfortunately, mean field theory can not tell reliably which 
of the possible competing instabilities wins. Therefore we are forced to look for finer 
methods, taking into account also the correlations which are neglected at the mean 
field level. The present paper is concerned with developing a systematic method that 
includes the dominant correlation effects neglected in the mean field theory. We require 
that the method satisfies the following criteria. 

(i) As emphasized by Anderson j^, the energy of a quantum state is dominantly 
determined by its short-range correlations. Therefore, an unbiased comparison of 
different symmetry breaking patterns requires that the same 'amount' of correlations 
is kept in all competing states. Obviously, mean field theory satisfies this criterion, 
neglecting the correlations in all states. 

(ii) Long ago it has been shown that, when treated beyond the mean field level, 
even purely repulsive systems can support superconductivity ■ We require that the 
general method reduces, in the case of superconductivity, to the recently developed 
perturbative scheme which enables a variational treatment of the Kohn-Luttinger 
superconductors 

Combining the above requirements, we are led to the search for a general 
lowest-order perturbative correction scheme to the mean field theory which treats 
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all symmetry breaking states on the same footing. This goal will be accomplished as 
follows. First we choose the type of instability we want to study. Then we split the 
Hamiltonian H of the interacting system into two parts. The first part contains the 
kinetic energy Hq — a ^k'T-k.o- and that part Hi of the interaction Hamiltonian 
which corresponds to a generalization of the reduced interaction Hamiltonian of 
the BCS theory 0]. For instance, in the case of superconductivity, Hi describes 
scattering of the Cooper pairs with zero total momentum, or in the case of magnetism, 
scattering of triplet particle-hole pairs with a given total momentum q. Denoting the 
rest of the interaction Hamiltonian as H2, the total Hamiltonian can be written as 
H = Hq + Hi + H2 ■ 

As shown in the case of superconductivity by Bogoliubov et al 5 , the mean 
field approach, when applied to the Hamiltonian Hq + Hi , is essentially exact in the 
thermodynamic limit. In order to preserve this exactness, the term H2, which describes 
scattering processes not included in the mean field theory, should be eliminated. 
Within our approach, this goal is achieved perturbatively. Namely, we construct a 
canonical transformation H = e^^He~^^ which eliminates scattering processes to first 
order in the interaction strength. This is achieved for H2+i[S, Hq] = and the effective 
Hamiltonian reads, to second order in the interaction strength, as H = Hq + Hi+ H2, 
where H2 — [iS,Hi + H2/2]. The term if 2, which has replaced the term H2 in the 
original Hamiltonian if, is typically even more complicated than H2. So what have we 
gained? The key observation of our approach is that at weak coupling, H2, being of 
second order in the coupling constant, is much smaller than H2 and therefore can be 
treated in the mean field approximation. Note that this yields a nontrivial correction 
to the original mean field theory, because H2 contains in general also terms of the 
type singled out into Hi. In order to better understand the last point, it is useful 
to observe that our approach is essentially variational, with a mean-field type ansatz 
for the wavefunction \^p) which minimizes the ground state expectation value of the 
energy, E = {ip\H\ip). Equivalcntly, one can write E = (V'l^^lV') with IV') = e'^'^IV') 
and interpret \ip) as a mean field state written in terms of quasiparticles. In this 
language it is not surprising that the effective interaction between the quasiparticles 
differs from the bare interaction. 

The plan of the paper is as follows. In Section 2 we construct the canonical 
transformation for the Hubbard model. We show explicitly that the effective 
Hamiltonian H depends on the instability channel. It is worth pointing out that 
this is a common feature of our method and of various versions of the rcnormalization 
group approach to correlated electrons P 171 151 lUl ITU] . 

In Section 3 we show how the method introduced in [3] fits the general scheme and 
we briefly review the main results of the method when applied to the superconducting 
instability. 

In Section 4 we develop the formalism for particle-hole instabilities with a finite 
total momentum of the particle-hole pairs. Our formulation allows for a simultaneous 
discussion of both the singlet and the triplet channels, the triplet channel being 
analyzed within an explicitly spin-rotation invariant formalism. 

In Section 5 we study particle-hole instabilities with a vanishing total momentum 
of the particle-hole pairs, called Pomeranchuk instabilities within Landau's Fermi 
liquid theory JT]. We take into account simultaneously the singlet and the triplet 
sectors, and also all point group symmetries, except for the identical representation 
in the singlet sector. For lattice models, this latter 'instability' is not associated with 
any symmetry breaking and simply corresponds to a change of the electron dispersion. 
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Finally, in Section 6 we present the results of explicit calculations for the t-t' 
Hubbard model on the square lattice at electron densities p corresponding to the Van 
Hove filling of the noninteracting band, and we specialize to p ~ 1 and small t'/t. 

2. Canonical transformation 

For the sake of simplicity, let us proceed with an explicit evaluation of H for the 
minimal microscopic model of interacting electrons, namely the Hubbard model. In 
that case we can write 

= XI 4T'^lT4iC2iAi234, 
{123} 

-f^2 = ^ cl^Cilcl^C2l (1 - A1234) , 
{123} 

where L is the number of lattice sites and the summation index {123} means 
momentum conservation ki + k2 = ks + k4. The cutoff function A1234 is always 
zero, except for the scattering processes which are to be singled out into Hi in which 
case A1234 = 1. For technical reasons, wc also take A1234 = 1 for processes which 
conserve the energy, i.e. if £1 + £2 = £3 + £4- This latter assumption will always be 
made in this paper and, for the sake of simplicity, it will not be mentioned explicitly 
in the formulas that follow. 

One verifies readily that the canonical transformation we look for has the following 
Hermitian generator, 

ill ^ 1 - A1234 t t 
^ {m} £1 + ^2 - £3 - £4 
It is worth pointing out that for 

A1234 = A2134 = Ai2435 (2) 

which will be always assumed from now on, the generator is explicitly spin rotation- 
invariant. A straightforward calculation shows that the eflfective scattering term H2 
is 

V- V- (1- Ai2:i4)(l + A,,,,,,A-) 



2L2 E E e^+e2-£3-£4 
{123} {a/37} 

{Sijcl-^Cal - <53a4|Ci|)4^C2J.4^C/34 + {S2Scl^C0l - S4l3cl^C2l)c!^^Ca'lcl-^Ci^ . (3) 

Note that H2 describes three-particle collisions and as such (together with 
terms generated from higher orders of perturbation theory) allows for more 
complicated patterns of symmetry breaking than simple particle-particle or particle- 
hole condensates. Here we do not consider such possibilities. 

In order to proceed it is useful to specialize at this point to a particular symmetry 
breaking channel. In this paper we will discuss three channels: (i) superconductivity 
with Cooper pairs with total momentum q = 0, (ii) particle- hole pairs with a finite 
total momentum Q, and (iii) particle- hole pairs with a vanishing total momentum 
q = 0. It is possible to apply our method also to superconductivity with Cooper pairs 
with a nonvanishing total momentum, but we have not done so. As will become clear 
later, in each symmetry breaking channel, there are singlet and triplet sectors in the 
spin space and various representations of the point group, therefore the number of 
different symmetry breaking states is enormous. 
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3. Superconducting channel 

This case has been discussed in detail in El El i so the discussion will be 
very brief. According to the BCS theory superconductivity can be viewed as an 
instability of the symmetric phase with respect to formation of Cooper pairs with zero 
total momentum. Therefore in the superconducting channel we choose 



^234 



for ki + ka = kg + k4 = 
otherwise 



Taking the expectation value of H in the superconducting state and introducing the 
angle-resolved pair field 6^/3 (k) = (c-kaCk/s), we have 

E = J2 ^k/k- - i E VipKpi^Mp) + i?FL, (4) 

ko" k.p 

where 

is that part of the interaction energy which is not associated with symmetry breaking. 
Note that we have chosen an opposite sign convention for V^^ with respect to 
El El El- In the present paper attractive interactions are chosen to be positive. 
For exponentially small order parameters, which we expect in the validity region of 
the present theory, we can neglect the change of -Bfl between the symmetric and 
symmetry broken state and this term will not be discussed any more. 
The pair scattering amplitude reads 

^^kp ^-U~ C/'xUlk + p, £p - Ek), 

where Xph(l''^) the real part of the particle-hole susceptibility 

/k - ./K^ 



Xph q, ^) = 7 V 



- EK-t-q -e-K-uj -iQ 

The imaginary part of Xph(q, cj) does not contribute to Vj^p, because the energy 
conserving processes are excluded from H2 ■ In the actual numerical implementation we 
neglect the feedback effects on the susceptibility and replace the occupation numbers 
/k by their values in the noninteracting system, This replacement is again well 
controlled for exponentially small order parameters. 

Following |15[ El we introduce a 2 x 2 matrix notation &k for the pair field with 
the matrix elements (b]s)af3 = ^a/3(k). Anticommutation of fcrmion operators then 
implies 6_k = ^^k where is a matrix transposed to X. Let us further define the 
gap matrix Ak — X]p I4;p&p with the same symmetry properies as &k and introduce 
the following parameterization 



Ak = 



-dl + idl, dl + 4 

di + tdi 



where a complex four-vector field d^ = (d^, dk) has been introduced which satisfies 
the conditions d'ij. = d^ and d_k = — dk- 

States for which the vector function qk = c^kd^ + ((i[^)*dk + idk x d^ vanishes 
identically are called unitary J31, since in that case AkAj^. is proportional to a unit 
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matrix. For unitary pairing states the theory simphfies considerably and the self- 
consistent equations for the gap function and the chemical potential /i, together with 
the expression for the ground state energy read as 



p 



2T 



(5) 



TV 



E 

k 



■ tanh 



^GS - - 



E 

k 



Sk 
(gk - ^k) 

2Ek 



2T 



where Ei^ = {£,1 + |Ak|2)i/2 ig the BCS quasiparticle energy, |AkP = J^l. 
the spectroscopic gap, is the bare electron dispersion, and ■Jk = £k — A*- 



\d^\^ is 



4. Density wave channel 

Next we turn to the discussion of instabilities of the symmetric state with respect to 
the formation of bound particle-hole pairs with a nonvanishing total momentum Q. 
The finite value of Q implies the presence of spatial modulations in the ground state, 
hence the name density wave channel. For the sake of simplicity, we specialize to the 
case when 2Q is an inverse lattice vector. In particular, this is relevant for the t-t' 
Hubbard model at the Van Hove density, when Q = {tt,tt); this latter case will be 
treated later as an explicit numerical example. 

The key quantities describing the symmetry broken phase are the angle- 
resolved particle-hole order-parameter fields dk = ^^(cj^^Cka) and dk = 
X]Q/3('^iiQ,'''c(/3Ck/3), describing pairs with spin zero and spin one, respectively. Here 
we have introduced the abbreviation k = k + Q and CTq/j are the Pauli matrices. Let 
us notice that = (d^)* and d^, = {d^)* . 

Next we ask the question about the effective interactions in the density wave 
channel. In order to single out into Hi those processes which scatter the particle-hole 
pairs with total momentum Q, we have to choose 

1 for ka — ki = Q or ks — k2 = Q 

otherwise 

Note that this choice leads to a spin rotation-invariant theory, since it satisfies the 
criterion Eq. [3 

Now let us calculate the expectation value of the Hamiltonian i/ in a density 
wave state |?Adw), E — {ipdw\H\ipdw)- After a tedious but straightforward calculation 
wc find the result 

E = Y, £k/k. -jY. [<^dUl + <"dk • dp] + ^FL, (6) 

kcr k.p 

in complete analogy with Eq. 01 The coefficients in front of the order parameter fields 
are the sought effective interactions. They read as 

T/kp^ =-[/+— [xpp(k + p, £k + £p) + Xpp(k + P, £k + £p)] 

- [x'phO^ -P,£k Sp) + Xph(k - P, ffk - £p)] , 
1/2 

Vkp^ - C/ - — [xpp(k + p, £k + £p) + Xpp(k + P, £k + £p)] , 



A 



1234 
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where XpplQiw) is the real part of the particle-particle susceptibility, 

Xpp(q,^^) 



1 1 - fO - 



L £k. + £q-K -UJ -iO' 

The imaginary part again does not enter due to the choice of A1234 and the actual 
occupation numbers have been replaced by their noninteracting values. Note that 
unlike in the superconducting case, the effective interactions are different in the singlet 
and triplet sector. 

Let us notice that both density wave interactions are real and have the following 
symmetries: Vkp = V^k — Vkp- analogy to the BCS theory, it is useful to introduce 
the gap functions A° = i"^ Ep ^kp'^'^p and Ak -^"^Ep^kp"4- O^e verifies 
easily that A? = (A^)* and Ajj — (Ak)*. Let us define the auxiliary quantities 
= (ck — £k)/2 and tJk = (sk + £k)/2 with the symmetry properties 6^ = —(5k and 
Wk = Wk, in terms of which the mean field Hamiltonian can be written in a compact 
form as 

H = ^(ck)t7VkCk + (^kAl^ + 4 • ^ 

k k 

where the prime restricts the summation to only within the magnetic zone and 
(ck)^ = (Ck|: cj^jj c^|i Cj^|) is a 4-component vector. Furthermore we have introduced 
a 4 X 4 matrix iVk — uj]^! + Mk, where I is a unit 4x4 matrix and 



M -(AO)*! -A* 

-A° 1 - Ak • -^kl 



with 1 denoting a unit 2x 2 matrix. Let us define the vector field Qk — AkAk + 
(Ak)*Ak — iAk X Ak- In analogy to the discussion of triplet superconductors, we will 
call states with Qk = as unitary. A major simplification is that for unitary states, 
is proportional to a unit matrix, = E^I with i?k = \/Sl. + |Akp, where 

we have defined the spectroscopic gap |Akp = E^=o l^kP- checks easily that 
the eigenvectors of Mk are E'k and —E^, both of them being doubly degenerate. In 
what follows we will develop an explicitly spin rotation-invariant mean field theory for 
unitary states, closely following the standard treatment of triplet superconductivity 

Let us assume that the unitary transformation Ck = C^k7k from the bare electrons 
described by (ck)^ to the new quasiparticles described by (7k)^ — ilkp^ti'^'t^'^ti) 

brings A^k to a diagonal form, iVk = ?7k^kC^k = diag(ii'k2, £^k2, £^ki, -Eki), where 
-E-ki = Wk + -Ek and £'k2 = '^-'k ^ -^k- Then the diagonalized Hamiltonian reads as 



[i?k27L7k. + EkilLlka] + E (^kA^ + 4 • A, 



kCT k 

and therefore (7k(7k)^) — 2^^ I + tanh ^iVk/2T^ j . Transforming back to the Ck 
operators, we thus have 



(ck(ck)t) = C/k(7k(7k)^)C/^ = 



(7) 
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Expanding the function tanhx into the Taylor series, making use of the binomial 



formula for powers of iVk 
find 



oJiJ. + Mk, and taking into account that AI^ 



tanh 



2T 



tanh 



tanh 



V 2T 



2Ek 



tanh 



Eki 
'2T 



tanh 



\ 2T 



Inserting this result to Eq. [7| and comparing the components of the matrices on both 
sides we find 



kl 



1 



dk =A? 



/kl 
/k2 



in 

Sk 

in 
Sk 



/i 



k2 



/k2 



1 



ill 

Sk 

in 
Sk 



/] 



kl 



■£;ki-£;k2' 



dl 



/k2 — /kl 



'^ii'i (cLck_o-> = 0. 

and /ki = f{Eki — fj,) where f{x) is the Fermi-Dirac distribution function. Note that 
/ki = /ki since i?ki = E^i. From here we finally find the self-consistent equations for 
the gap matrix and for the chemical potential, as well as the corresponding ground 
state energy: 

1 ,,=.^„, /p2 ^ /pi 



-C/kl — ^k2 

In the above equations we have used the notation rikcr — (cjj.^Cko-) 



Ak 



^ T/sdw /\ 



1 



A^ =t'E<"(A°) 



p2 

/p2 — /pi 



L ^ 
p 



pi 



E, 



p2 



(8) 



(9) 



k 



[/kl + /k2] 



Skl/kl+i;k2/k2 + |Ak| 



/k2 - /kl 
_Ekl — £^k2 



Note that the summations in these equations run over the full Brillouin zone. It should 
be stressed that the simple BCS-like form of Eas. I8I9I obtains only for unitary states 
with Qk = 0. 

In what follows, we do not allow for the simultaneous presence of charge and spin 
density wave order. Nevertheless, even with this simplification the number of possible 
phases turns out to be quite large. For instance, on a square lattice there are ten 
different symmetry breaking patterns in the singlet sector, since A^ may transform 
according to one of the five irreducible representations of the point group and there is 
an additional double degeneracy associated with the parity of the gap function under 
the translation in momentum space, k ^ k. Namely, if we decompose the order 
parameter to its real real and imaginary parts, 
be written as 



A^ 



2;k + ij/k, the gap equation can 



Xk 

2/k 



1 T/cdw 
P 



fp2 - /pi 

Epi — E, 



p2 



As observed by Hankevych and Wegner |1U| . the components with the symmetry 
properties x^^ = Xk and ^k = — j/k are described by coupling constants with opposite 
signs. 
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As regards the triplet sector, based on the analogy of Eq. |S1 with the triplet 
superconductor case ^1 we expect that for degenerate representations of the 
point group, energy is minimized by Balian-Werthamer-like states involving complex 
patterns of Ak- For square lattices only the p representation is degenerate and since 
that symmetry sector is not important for the model studied in Section 6, in this paper 
we discuss only states of the type Ak = (a;k + iyk)n^ where n is a fixed direction in 
spin space and Xk, j/k are real functions with the symmetry properties = a;k and 
2/k — — 2/k- Repeating the argument for the charge density waves, we find again that 
the components Xk and j/k are described by coupling constants with opposite signs. 



5. Landau channel 



In this Section we consider instabilities of the high temperature state with respect 
to the formation of bound particle-hole pairs with total momentum q = 0. In this 
so-called Landau channel we therefore choose 

1 for ka — ki or — k2 
otherwise 

Let us introduce the angle resolved singlet and triplet particle-hole fields, rik — 
2~i X]a(cLcko-) and dy, = '^'^J^api^la^afiCkp), respectively. In terms of rik and 
dk, the expectation values of Hi and H2 read as 



A1234 ~ 



(Hi) = ^ ['^k?^p - dk- dp 
kp 



/T~r \ Sr^ 1-711-712 

{H2) =J2 2^ 



713714 — d^ ■ di 



A crucial difference between the singlet and triplet fields is as follows. The triplet 
fields dk vanish in the symmetric high-temperature phase and we expect that, in the 
weak coupling limit, the development of a finite dk does not substantially change 
the electron distribution function. Therefore we can construct the triplet Landau 
interaction function as a second derivative of the interaction energy (Hi + H2) with 
respect to dk- The derivatives can be evaluated for undeformed Fermi surfaces, leading 
to the result 

^kp ~ C/2xpp(p -f k, Ep + £k). 
Contrary to the triplet case, the singlet field Uk does not vanish even in the 
noninteracting system. If we expand the interaction energy {Hi + H2) with respect to 
deviations of d'^ from its value in the noninteracting system, besides quadratic terms 
there appear also linear terms in the deviation. Within the Landau Fermi liquid 
theory, this means that the quasiparticle dispersion relation is modified with respect 
to the bare spectrum. However, since in the discussion of superconductivity and of the 
density waves we have not taken into account the renormalization of the spectrum, 
in order to keep the same level of approximation, we will consider only those singlet 
sector Landau instabilities which occur in the non s-wave channel, i.e. we assume 
that 71k = /k + c^k where is the noninteracting distribution and d'^ transforms with 
respect to a nontrivial representation of the point group. With this restriction the 
terms linear in d'^ vanish and we can write 

E = Y^ ek/k. -^Y. + ■ dp] + Efl, (10) 

kfT k,p 
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where 



^kp - [x^plP + k, ep + £k) - 2x;h(P - k, ep - ek)] 



Note that both in the singlet and in the triplet sectors, the interaction matrix is 
real and symmetric, Vkp = Vpk- After introducing the real gap functions Aj^ = 
J2p l^kpC^p E^nd Ak = L^^ X]p ^u.p'^P' mean field Hamiltonian can be written 

as 

H = 5](ck)t [(£k - AS^)1 - Ak • a] ck + [d^Al + dk • A 



where (ck)^ — (Ck|,cj^|)- Consider a unitary transformation Ck ~ C^k7k to the 
quasiparticle operators (7k )^ = (7k-i7k+) with 



C/k 



1 



2IA1 



Ak| + ^ 

\/|Akl+AJ 



where we have introduced |Akp = ELi(^k)^ 
Hamiltonian is diagonal in the new basis. 



^kl-A^ 

Vl^kl-Aj 

One checks readily that the 



H = 



J2 [-Sk-7L7k- + -Bk+7i^+7k+ + 4K + d^. ■ Ak 



with energy eigenvalues i?k+ = £k — A[^ + |Ak| and E^k- = Ck — A[^ — |Ak|. Repeating 
the argument leading to Eq. Owe find 

1 - /k- 



(ck(ck)^) = f/k 







l-/k 



Making use of the formula J/kC^C/j^ = (T- Ak/| Ak|, the right hand side can be computed 
explicitly. Comparing the matrix elements of the left and right hand sides we find 



Ai 



/k- - /k+ 



4 



1 



[/k- + /k+] - /k 



From here follow the self-consistent equations for the gap functions and the chemical 
potential, as well as the corresponding ground state energy: 

Ak - y2^^^kpAp — 



p+ 



^p- 



A^ 



N 



(/p- + /p+) - fp 



(11) 
(12) 



Egs = 





Al!\ 








+ ^)/k-^ 


- (i^k+ 





/k+ + 



IAJ2 



/k 



-Ek+ — i?k- 



We stress that Eas- llllf^ applv only if singlet channel Pomeranchuk 'instabilities' are 
not allowed. In what follows, we specialize to pure singlet or pure triplet instabilities. 
Note that, as in Section 4, for degenerate representations nontrivial states in the triplet 
channel, in which Ak does not preserve the same direction in spin space as one varies 
k, are likely. 
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6. Application: t-t' Hubbard model at the Van Hove density 

In the rest of this paper, we will study the phase diagram of the square lattice t-t' 
Hubbard model in the plane with the ordinates t' /t and electron filling, p. In order to 
be able to discuss also the density wave states, we restrict ourselves to the so-called 
Van Hove line, i.e. to those combinations oit' /t and p, which lead to a Fermi surface of 
the noninteracting system which crosses the saddle points of the dispersion at (±7r, 0) 
and (0, ±7r). With this choice the most singular density wave susceptibility is expected 
at Q = (tTjTt), satisfying the criterion that 2Q is an inverse lattice vector. 

All results reported in this Section were obtained on L x L lattices with L — 128 
sites and periodic boundary conditions. The susceptibilities entering the effective 
interactions have been written as convolutions and calculated using the Fast Fourier 
Transform algorithm, as explained for the particle- hole case in JJ^. In the time 
direction, we have sampled the interval (0, Tmax) with N points. Thus the energy 
precision is Aw = 27r/Tmax and the maximal frequency is Wmax = ^wN/Tmax- In order 
that Auj is comparable to the energy precision in k-space, we require Aoj = 8t/L. 
We also require Wmax = 64t in order to faithfully describe the high energy processes. 
This leads us to the choice N — 8L ^ 1024 time points and Tmax = 7rL/(4i) = 327r/t. 
In order to reduce the numerical error, the time evolution was damped according to 
exp(— Ft) with F — Auj/2 = t/32. The accuracy of the FFT algorithm was checked 
by direct calculation of the susceptibilities. The self-consistent equations were solved 
by the damped iterative method. The condensation energy i^cond = Enoim — Eqs was 
calculated as the difference between the energy of the symmetry broken state Eqq and 
the normal state energy i?norm = 2X]k^k/k- Note that stable states are described by 
a positive i^cond- 

In what follows the symbols SC, CDW, SDW, SL, and TL stand for the 
superconducting, charge density wave, spin density wave, singlet Landau channel, and 
triplet Landau channel phases, respectively. Spatial representations of the point group 
of the square are denoted s, d, d^y, and g (one-dimensional, even representations), and 
p (two-dimensional, odd representation). Note that we have not introduced separate 
symbols for singlet and triplet superconductors, since for superconductors the spin 
sector is completely defined by the parity of the spatial representation. The density 
wave instabilities are characterized by an additional quantum number |10| . namely 
the parity of the order parameter under translation in momentum space by Q. Order 
parameters satisfying Aj^. = Ak and A^ = — Ak are distinguished by the suffix 1 and 
2, respectively. 

In Fig.Qwc plot the condensation energy as a function of U for t' /t = and at half 
filling, p = 1.0. As expected, at weak coupling the leading instability is towards an s- 
SDWl (antiferromagnetic) instability. Around U ^ it the condensation energy of the 
antiferromagnetic state acquires a maximum. At larger coupling constants different 
symmetry breaking patterns dominate. Since the Hubbard model at half filling and 
t' /t = is expected to order antiferromagnetically for all coupling constants J^, from 
Fig. n we estimate that our calculations are qualitatively correct up to C/ w ?>t. 

It should be noted that the overall shape of Fig.^ and in particular the prediction 
of the relative stability of the phases, is in qualitative agreement with the calculation 
of the transition temperatures by the flow equation method by Hankevych and Wegner 
|1()| . see their Fig.l. The only qualitative difference with respect to jTO] regards the 
relative stability of the d-LS (d-wave Pomeranchuk instability ^11201) and p-CDW2 
(band splitting ^U]) phases: according to our calculation, the p-CDW2 phase is more 



Phase diagrams of correlated electrons: systematic corrections 



11 




Figure 1. Condensation energy as a function of U at tlie Van Hove filling, 
t'/t = 0, p = 1.0, L = 128. The d-SC state is degenerate with d-SDW2 and 
d-CDW2 states. 




Figure 2. Gap function in the antiferromagnetic (s-SDWl) phase. The 
parameters are the same as in Fig. Q and U = 3t. 

stable than the d-LS phase. Hankevych and Wegner have also observed that for 
t'/t = and p — 1.0, there is an additional symmetry which guarantees that the d-SC 
state is degenerate with the (i-SDW2 (triplet flux phase j21j) and c?-CDW2 (singlet 
flux phase |^) states. This degeneracy can be proven also within our formalism and 
is nicely satisfled by the data, thus providing a nontrivial check of the numerics. 

Before proceeding it is worth pointing out that the gap function in the simple 
antiferromagnetic (s-SDWl) phase is by no means featureless, as it would be in a 
mean fleld theory. Figure |5] shows that Ak is suppressed along the underlying square 
Fermi surface, the suppression being largest close to the saddle points at (±7r, 0) and 
(0, ±7r). This suppression is caused by the repulsive term in the effective interaction 
Vjlp*. Figurenshows that the repulsive term starts to dominate at [/ > 3t. A proper 
discussion of the spectroscopic implications of Fig. |21 would require including also the 
change of the normal state dispersion (s-SL channel), which we have not attempted. 

Summarizing the results at half filling, our data indicate that already at 
U ~ 3t where our approach should still apply, there exist nonnegligible tendencies 
towards ordering in several nontrivial symmetry breaking states, many of which have 
been extensively discussed in the context of the physics of the cuprates: d-wave 
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Figure 3. Condensation energy as a function of U at tlie Van Hove filling, 
t'/t = 0.05, p = 0.959, L = 128. 




Figure 4. Condensation energy as a function of U at the Van Hove filling, 
t'/t = 0.10, p = 0.918, L = 128. 




Figure 5. Condensation energy as a function of U at the Van Hove filling, 
t'/t = 0.15, p = 0.875, L = 128. 



superconductivity, d-wave Pomeranchuk instability, and the flux phases (in the singlet 
case also called d-density wave phase [221 )■ what follows we discuss the evolution 
of these tendencies as we move along the Van Hove line towards smaller electron 
densities. 

Inspection of Figs. and I3I5I shows that the antiferromagnetic (s-SDWl) 
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0.865 0.87 0.875 0.88 0.885 0.89 

rho 

Figure 6. Condensation energy as a function of p in the vicinity of the Van 
Hove fiUing, t' /t = 0.15, U = 3t, L = 128. 

condensation energy rapidly diminishes with decreasing fiUing. This opens up the 
interesting possibiHty that a nontrivial phase is stabihzed at the Van Hove hne. The 
leading subdominant instability at half filling is towards the band splitting phase 
(p-CDW2). Figures n El and 21 show that this phase quickly looses stability and 
never becomes dominant in the region of applicability of our method. Moreover, 
the degeneracy between the d-SC and the flux phases is lifted away from half filling: 
i;cond(dSC) > £;cond(dSDW2) > £;cond(dCDW2). 

Thus we are led to a study of the competition between the antiferromagnet, 
the d-wave Pomeranchuk phase, and the c?-wave superconductor. We find that, 
along the singular Van Hove line, the Pomeranchuk instability is stronger than the 
superconducting instability in the whole region we have studied, i.e. up to t' /t = 0.2. 
Moreover, if we fix J7 = 3i and consider the competition with the antiferromagnetic 
phase, we find that the d-wave Pomeranchuk phase becomes stabilized for t'/t > 
(t'/t),^ 0.10. 

At the mean field level, the system at the Van Hove density is unstable towards 
antiferromagnetism for infinitesimal coupling. This raises the following interesting 
question: Is it so that at sufficiently weak coupling, when the mean field theory should 
apply, the system always has to order antiferromagnetically? Figures □ andHagree 
with this hypothesis. At smaller electron fillings we cannot confirm it, since for our 
lattice size we cannot test reliably states with small condensation energies. 

Our results indicate that, at the Van Hove filling, the particle-hole instabilities are 
stronger than the superconducting ones. In the rest of this Section, we wish to study 
the relative stability of phases away from the Van Hove filling. First we notice that, 
in regions of the t'/t vs. p plane lying close to the Van Hove line, the density wave 
instabilities are expected to occur at wavevectors slightly different from Q = (7r,7r), 
and the formalism developed in Section 4 does not apply. Therefore, in order to work 
in a region where the density wave instabilities are safely negligible, we have chosen 
to study the competition between the d-wave Pomeranchuk phase and the d-wave 
superconductor in the vicinity of the Van Hove point at t'/t = 0.15 for U = 3t. The 
results are shown in Fig. El As expected, sufficiently far away from the Van Hove 
filling (located at p = 0.875), the superconducting phase is stabilized. An unexpected 
feature of the data in Fig. is that the d-wave Pomeranchuk phase is stable in an 
extremely narrow region around the Van Hove line even at moderate coupling U ~ 3t. 
This might explain why some of the renormalization group studies observed the d-LS 
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phase JHI while others did not [HI- 
7. Conclusions 

In this work we have introduced a simple scheme which enables us to construct 
effective Hamiltonians for correlated electron systems in a controlled perturbative way. 
The method has been applied to the Hubbard model, but more complicated models 
can be also studied. As it stands, our scheme allows us to treat superconductivity, 
Pomeranchuk instabilities (except for the singlet instability transforming according 
to the trivial representation of the point group), and density wave instabilities with 
ordering wavevector Q (subject to the restriction that 2Q is an inverse lattice vector). 
Both the singlet and the triplet sectors within each channel have been treated and 
it has been observed that, for degenerate representations, the phenomenology of the 
triplet particle-hole channels might be as rich as that of the triplet superconductors 

An effective implementation of the method, making use of the Fast Fourier 
Transforms, has been developed. Lattices as large as 128 x 128 are tractable on PCs. 
The method has been applied to the phase diagram of the t-t' Hubbard model on the 
square lattice close to the Van Hove filling, for t' /t in the interval between and 0.2. 
Our results are similar to those obtained by the flow equation method ^ni^ at the Van 
Hove line, particle-hole instabilities are found to dominate. The antiferromagnetic 
and the c?-wave Pomeranchuk phases are stable for t' /t < {t' /t)c and t' /t > (t' /t)c, 
respectively. The critical value of {t'/t)c depends on the value of the interaction 
strength U. The particle-hole instabilities are confined to an extremely narrow stripe 
of the t' /t vs. p plane around the Van Hove line even at moderate coupling strength 
U = 3t, leaving most of the phase space to be dominated by the Kohn-Luttinger effect 

m 

Straightforward future improvements of this work might include the study of 
larger systems and a finite size scaling analysis of the results. The theory can be 
easily extended to finite temperatures. It is also possible to take into account the 
self-energy effects, or, in other words, the s-LS channel, and also the change of Epi^ in 
the symmetry breaking phases. Furthermore, the method can be applied to density 
wave states with more general wavevectors or to states with simultaneous presence of 
symmetry breaking in different channels or symmetry sectors. 

Turning to more speculative issues, an intriguing generalization would be to 
consider condensates of more than two particles or holes. A bosonic condensate would 
require a quadruplet of fields and an inspection of Eq.O suggests that in order to find 
the effective interaction for such processes, the perturbation expansion would have to 
be continued to higher orders. Finite lifetime effects should be studied as well. This 
might require a reorganization of perturbation theory in order to take into account 
the fact that the energy conserving processes have been singled out into the tractable 
term Hi. 

A serious limitation of our scheme is its inability to treat states without obvious 
order parameters, such as the Mott insulators. In such cases the dynamical mean field 
theory (for a review, see |2S1) seems to be the method of choice. On the other hand, 
for problems requiring a high resolution in momentum space, our method might be 
preferable. 
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